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A correlated random walk approach to diffusion is applied to the disordered non-overlapping 
Lorentz gas. By invoking the Lu-Torquato theory for chord-length distributions in random media 
[J. Chem. Phys. 98, 6472 (1993)], an analytic expression for the diffusion constant in arbitrary 
number of dimensions d is obtained. The result corresponds to an Enskog-like correction to the 
Boltzmann prediction, being exact in the dilute limit, and better or nearly exact in comparison to 
renormalized kinetic theory predictions for all allowed densities in d = 2, 3. Extensive numerical 
simulations were also performed to elucidate the role of the approximations involved. 
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I. INTRODUCTION 

The problem of computing the transport coefficients 
of coUisional models from first-principles can be ap- 
proached with a variety of methods, each with differing 
scope and complexity. At the most elementary level, ki- 
netic arguments based on the concept of mean- free-path 
yield estimates with the correct dependence on the rele- 
vant parameters, but as one might expect miss the cor- 
rect prefactors even at low densities |T|. Progress can 
be made using the Boltzmann equation and the asso- 
ciated Chapman-Enskog method, which give the exact 
low-density limit, though at the expense of substantially 
greater effort |2]. At higher densities, renormalized ki- 
netic theory accounts for the peculiar divergences one 
faces when expanding the coefficients in powers of den- 
sity ^, though here too the improvement comes at the 
cost of increased complexity in the calculation. 

For the case of the periodic Lorentz gas - i.e. the 
simple model of a point particle colliding with a fixed 
array of spherical scatterers - an alternate approach of 
great simplicity that bypasses the above kinetic route 
has been proposed by Machta and Zwanzig [4 . These 
authors recognize the fact that for small enough spac- 
ings between the scatterers (i.e., in the opposite limit of 
the dilute approximation), the diffusing particle spends 
most of its time trapped inside effective "cages" before 
finding its way into an adjacent trapping region. The 
particle motion can thus be approximated as a random 
walk on a lattice, where the rate of transition is given by 
the cage escape rate, and consequently a simple analytic 
formula for the diffusion constant D can be derived. The 
method has been subsequently extended and refined (see 
e.g. m m), but its success remains contingent upon 
the prevalence of such trapping events. 

The approach adopted in this paper — though concep- 
tually different from that of Machta and Zwanzig — is an 
attempt to arrive at a result of similar simplicity and 
accuracy for the disordered case, without relying on the 



FIG. 1: (Color) A typical trajectory of the non-overlapping 
Lorentz gas in d = 2 at packing fraction (p = 0.4. The point 
particle propagates freely at constant velocity v until it suffers 
specular collisions with the scatterers, each of radius R. Inset: 
Representative segment of the trajectory trapped in an effec- 
tive "cage," an event that is ultimately responsible for corre- 
lations in the displacements Aa; between collisions. In higher 
number of dimensions or at lower densities, where greater 
voids exist between neighboring scatterers, such events are 
effectively suppressed. 



•Electronic address: |adiba@mail.nih.gov| 



existence of kinetic traps. The method is a simple appli- 
cation of correlated random walks (CRW), whose main 
result is the ability to express D very simply in terms of 
the correlation and moments of free-paths [8]. For the 
present problem, these can be easily computed by invok- 
ing the scattering properties of point particles by hyper- 
spheres (Eq. (|4|) and the distribution of chord- lenghts in 
random media 0, respectively. A closely related apph- 
cation of CRWs has been used in the study of Knudsen 
diffusion in porous media \iO\ 1111 . In the following, the 
correlated random walk method for the computation of 
the diffusion constant will be briefly reviewed and applied 
to the d-dimensional Lorentz gas. All calculations will be 
reported for the disordered non-overlapping case, where 
the scatterers are distributed like a liquid in equilibrium. 



II. THEORY 

To begin the derivation, note that the net displacement 
of a particle starting from the origin consists of discrete 
steps Axi between successive collisions i — 1 and i (see 
Fig.ll]), so that x(nc) = X^T^i ^'^ii where Uc is the num- 
ber of collisions. Thus, using the fact that Uc is related 
to the continuous time t through t — X]"^! ^^' ~ ""c (At) 
(the last equality being valid in the large Uc limit) and 
that Ai = Ax/v where v = |v| is a constant, routine 
manipulations of the mean square displacement (x^(nc)) 
give an exact formula for D expressed solely in terms of 
the statistics of the collision vectors, namely 



D 



2d(Aa 



(Ax2)+2^(Axo-Ax,) 



1=1 



(1) 



where Axq is the displacement of the particle before its 
first collision. This expression can be seen as the discrete 
analog of the so-called Green-Kubo formula for the dif- 
fusion constant in terms of the continuous velocity corre- 
lation function (curiously, the latter was already derived 
by Taylor himself decades before Green and Kubo [S], 
suggesting the more appropriate nomenclature Taylor- 
Green-Kubo for that result [7]). 

The assumptions specific to the Lorentz gas that will 
be made (and tested) in the remainder of the paper can 
be cast directly in terms of the statistical properties of 
the collision vectors Ax = Axv, where v is the unit 
velocity vector. To wit, it will be explicitly assumed 
that (a) their magnitudes Ax are uncorrelated, and that 
(b) their orientations v change in a Markovian fashion. 
The last assumption is the simplest condition consistent 
with the dynamics of the gas particles, as the underlying 
specular collision rules manifestly introduce angular cor- 
relations between successive velocity vectors, and hence 
these cannot be assumed to be completely uncorrelated. 
As such assumptions are essentially equivalent to Boltz- 
mann's molecular chaos hypothesis, they are satisfied for 
a random arrangement of scatterers at sufficiently low 
densities, and thus the ensuing predictions for D should 
be exact in this limit. 

With assumption (a), the infinite sum in Eq. (fTl) sim- 
plifies to (Ax) J2hLi {^0 ' ^i) ^ ^^'^ ^^ accordance with 
(b), we can invoke the Chapman-Kolmogorov relation to 
express the conditional probability P(vi|vo) in terms of 
P(vj|v;_i) and hence obtain the recursion relation 

(2) 



(Vo • V;) 



(Vq •V;_i; 



that 



where {cos 9) = (v/ • v;_i) (this result assumes 
P(v;|vi_i) has azimuthal symmetry about Vi_i, as in 
the case of scattering by hyperspheres). This relation 
can be iterated to yield the simple result [19] 



(Vo ■ V;) 



(3) 



Now, from the specular collision rules of the gas, a direct 
geometric computation assuming a uniform flow of inci- 
dent particles towards a d-dimensional spherical scatterer 




FIG. 2: (Color) Discrete velocity correlation function of the 
Lorentz gas at packing fraction (ji = 0.1 in d = 2 and d = 3, 
according to simulation results (main panel, symbols) and 
the Markovian prediction given by Eqs. (|3| and Q (main 
panel, lines). Inset: Relative error between simulation and 
the Markovian prediction for the integrated autocorrelation 
function Y.'iZo ('^o • v;). 



gives the following expression for the average scattering 
cosine: 



d-3 



cos W = 



(4) 



a result that can be easily checked against the usual dif- 
ferential cross-sections of hard disks and hard spheres 
(compare with the analogous result for the case of porous 
media in d = 3 [10l[TT]). In Fig. [2] the above predictions 
for the discrete velocity correlation function (Eqs. (pi) and 
(E|) are compared with simulation results, illustrating in 
particular the robustness of the Markovian approxima- 
tion even at high densities. 

Use of the foregoing in Eq. (IT|) reduces the dependence 
of D on the series of collision vectors to their first two 
marginal moments, namely 



D = 



v{Aa 



' (Ax^) 
2{Axf 



- 1 



(5) 



This result is a straightforward consequence of assump- 
tion (b) in correlated random walks (analogous to what 
was done originally in one dimension [S]), where the 
main additional ingredient is the explicit computation 
of the degree of correlation between adjacent v, i.e. 
(v; -v;-!) = (cosO), givcu by Eq. ^ for the present 
problem. Generalization to non-spherical obstacles is 
easily performed, provided one knows their scattering 
properties so that a formula similar to Eq. Q can be 
obtained. 
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FIG. 3: (Color) Distribution of collision distances p{Ax/R) 
in the Lorentz gas for d — 2 and d = 3. The distributions 
are arbitrarily normalized for clarity of presentation. The 
symbols are simulation results at the indicated packing frac- 
tions, while the lines are exponential distributions with the 
exponents A predicted by the Lu-Torquato theory for chord- 
lengths (Eq. Q). 



A. Low-density behavior 

To check the above prediction for D against known 
results in the dilute (Boltzmann) limit, first note that 
in this case the Ace are exactly exponentially distributed 
jT], and hence the term in parenthesis in Eq. (pi vanishes. 
The first moment of this distribution is given by l/{crp), 
where a is the total cross section of the obstacles and 
p = N/V is their number density. For d = 2 and d = 3, 
this yields, respectively. 
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where (p is the packing fraction of the obstacles, and the 
subscript emphasizes that these results are only valid in 
the Boltzmann limit. These are precisely the same pre- 
dictions obtained via more sophisticated methods based 
on, e.g., the Chapman-Enskog method or Zwanzig's op- 
erator expansion scheme ^121 1131 114] . In accordance with 



FIG. 4: (Color) The diffusion constant of the Lorentz gas 
in d = 2 (top) and d = 3 (bottom) according to the vari- 
ous theories mentioned in the paper. Legends are the same 
for both figures. The elementary kinetics predictions are 
D/{vR) = (7r/4)(/)-^ for d = 2 and D/{vR) = (4/9)0~^ for 
d = 3, the latter being identical to the Boltzmann predic- 
tion and hence not shown. The long-dashed lines are the 
Boltzmann results given by Eq. H, while the red curves are 
from renormalized kinetic theorv [T2l 113] . The analytic re- 
sults of this work (blue curves, closer to simulation results) 
are given by Eq. ([9|, and are visually indistinguishable from 
what one obtains by using Eq. ([5| with the required moments 
measured from the simulation itself (not shown). The arrows 
indicate the maximum packing fractions 0niax ~ 0.907 and 
9!>max ~ 0.740 for d = 2, 3, respectively j^. 



the general rule [Ij, such predictions differ from those of 
elementary kinetic theory — namely, Dj^ = v (Ax) / d — 
by factors of order unity. (An accidental exception in 
the present problem is the d = 3 case, due to the ef- 
fectively randomized nature of its scattering, as can be 
seen by the vanishing of Eq. (|4|). In terms of the present 
framework, it can be easily seen that the source of error 
in elementary kinetic theory lies in its complete neglect 
of correlations in v; indeed, had we made this assump- 
tion in our derivation, our expression for D would have 
been coincident with Dk- 
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FIG. 5: (Color) Correlation times of the collision distances 
Ax for d = 2,3 at various packing fractions, measured by 
summing tiie numerical results for the autocorrelation func- 
tion of Ax over all possible lags (/ = l,2,...,oo). Inset: The 
autocorrelation function itself for three densities in d = 2. 



B. High-density extension 

To go beyond the Boltzmann limit, a theory for the 
distribution of Aa; at high densities is necessary. Here 
we invoke the Lu-Torquato theory for chord-length dis- 
tributions in random media [T51 [1^] . This distribution is 
defined as the probabiHty of finding an unobstructed line 
segment of length A a; connecting the surfaces of any two 
obstacles [5]. The Lu-Torquato approach is conceptually 
simple and elegant, and uses the scaled-particle method 
of liquid theory to relate chord-length distributions to the 
problem of finding a thin cylindrical cavity in the liquid. 
The central prediction for non-overlapping spheres is that 
the distribution is exponential, viz. p(Aa;) = Ae" 
with exponent given by [9] 



-AAk 
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where ujd = 7r''/^/r(l + d/2) is the volume of a d- 
dimensional unit sphere. A comparison of these predic- 
tions with the present Lorentz gas simulations is shown in 
Fig.[3J With the exception of the salient non-exponential 
behavior for (/) = 0.8in(i = 2 (due to excessive trapping), 
the theory is in excellent agreement with the numerical 
results. 



Finally, with Eqs. (jsl and (l7| we arrive at our result 
for the diffusion constant of the d-dimensional Lorentz 
gas: 



D _ d+l uod 1 
vR Ad uJd-i 



(8) 



which for the particular cases of d = 2 and d = 3, reads, 
respectively, 

D ill 1-4) D 41-0 

Note that, incidentally, this result can be obtained 
through an ad hoc excluded- volume (i.e. Enskog-like |17j ) 
correction to the Boltzmann collision frequency, which 
ultimately adds an extra factor of 1 — to Eq. (l6| . 

The above predictions are compared with simulation 
as well as with the renormalized kinetic theory of Weij- 
land and van Leeuwen |T^ |T5] in Fig. [4] As immediately 
apparent, the theory is in excellent agreement with sim- 
ulation for d = 3, while the agreement for d = 2 is less 
satisfactory; this can be traced back to the sensitivity of 
assumption (a) to the dimensionality of the problem (see 
Fig.l5|, as already anticipated in Fig. fl] 



III. CONCLUSIONS 

To summarize, in this paper a simple analytic expres- 
sion for the diffusion constant of the disordered non- 
overlapping Lorentz gas in arbitrary number of dimen- 
sions was derived (Eq. (Is])). This result was obtained 
through a simple application of correlated random walks, 
combining the analytic predictions of Lu and Torquato 
(Eq. (|7|) with the exact scattering properties of point 
particles by hard hyperspheres (Eq. (B). The simplic- 
ity and accuracy of the ensuing theory in comparison to 
previous renormalized kinetic methods |12l 113] are re- 
markable (cf. Fig.l4|, and in this regard the goal set out 
in the introduction was achieved [20]. It is hoped that 
these encouraging results will foster the application of 
similar ideas to more general models of diffusion. 
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